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<D . Abstract 

We consider properties of the structure functions of inclusive heavy meson decays B — > 
X c and treat the c quark mass as a free parameter. We show that in two extreme cases 
^ of heavy and light c quark the structure functions of heavy-heavy and heavy-light transi- 

tions are given by a Fourier transform of the matrix elements of Wilson lines containing a 
time-like and a light-like segment, correspondingly. Using the renormalization properties 
of Wilson lines we find the dependence of the structure functions on the factorization scale, 
the structure function of heavy-heavy transition is renormalized multiplicatively while that 
of heavy-light transition obeys the GLAP-type evolution equation. We propose a gener- 
alization of the sum rules for the moments of the structure functions (Bjorken, Voloshin, 
and the "third" sum rules) with a soft exponential factorization cut-off, which correctly 
incorporates both perturbative and nonperturbative effects. We analyze nonperturbative 
corrections by first considering infrared renormalon contributions to the Wilson lines. Un- 
certainties induced by the leading renormalon pole at u = ^ are exactly cancelled by the 
similar uncertainty in the heavy quark pole mass. The leading nonperturbative corrections 
associated with the next renormalon at u = 1 are parameterized by the matrix element /j, 2 
which is proportional to the heavy quark kinetic energy. 
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1. Introduction 

Recently, the interest to the inclusive semileptonic decays of the B meson was renewed after it 
was proposed [1], [| to apply the operator product expansion (OPE) in powers of a heavy quark 
mass to the analysis of the heavy meson decays. Being combined together, the OPE and heavy 
quark effective theory allow us to expand the total decay rates of the B meson in powers of A/m^, 
with A the QCD scale and nib the b quark mass, and parameterize nonperturbative corrections 
by hadronic matrix elements |3| . The analysis of the differential distributions of the B meson in 
the framework of the heavy quark expansion turned out to be more complicated [|J. Let us take 
as an example the inclusive decay B — > W + X C with W —>■ W and consider the decay distribution 
dT-B^w+Xc/ dW 2 dWo with respect to the invariant mass W 2 and the energy Wo of the lepton pair 
in the rest frame of the B meson. The process has two large momentum scales, the heavy quark 
mass, mi,, and the invariant mass, mx, of the final state 

m 2 x = (m B v - W) 2 = m 2 B - 2m B W Q + W 2 

with v M the velocity of the B meson. The OPE of the differential distribution j|] is an expansion 
in powers of the b quark mass, A/m b , and virtuality of the outgoing c quark, Am h /{m 2 x — m 2 ), 
provided that both parameters are small. Although this expansion is well defined for large 
invariant masses mx, it becomes divergent in the region where 

m 2 x — m 2 c < Ami, ■ (1) 

In this region the OPE fails and the differential distributions in the decay B — > W + X c are 
parameterized in the leading A/mf, order by the universal nonperturbative structure function 
of the B meson B. This function takes into account large perturbative and nonperturbative 
corrections to the differential rate of the decay in the region (H). The structure function cannot 
be calculated at present out of QCD, but it turns out to be possible to relate its integral 
characteristics to a few phenomenological parameters. The corresponding relations are known as 
Bjorken || 0, Voloshin || [|, and the "third" flO| sum rules. They express first three moments 



of the structure function in terms of the fundamental parameters of the HQET. It was realized, 
however, |7|, [|, ||, that these parameters become functions of the normalization scale after we 
take into account perturbative corrections and the sum rules should be properly modified to 
incorporate perturbative effects. It remains unclear how to renormalize consistently the sum 
rules for the structure function and this is one of the questions we consider in the present paper. 
The properties of the structure function crucially depend || on the ratio of the quark masses 
m c /mb. We study below the structure function of the B meson and treat the c quark mass as a 
free parameter. 

It is convenient to define the angle, ■&, in Minkowski space-time between 4-velocities v = 
Pb/i^b and v' = Px/f^x of the B meson and X c , respectively, as ch?9 = (v ■ v'), and perform 
analysis using the variables $ and mx, 

W 2 = m x + m 2 B — 2m x r m B ch $ , W = m B — m x ch ■& . 

Then, the phase space for m x and $ is 

mo < mx < m B , < $ < In , 

m x 



where the maximum value of *& is achieved in the exclusive D meson production, mx = vtid. In 
this paper we study the B meson decay in the subregion of the phase space defined in ([[]). 

Considering all possible values of m c in ([I]), we distinguish two special cases of "heavy" and 
"light" c quark mass, m 2 ^> Am b and m 2 <C Am b , respectively. The real masses m b , m c lie 
somewhere between these two extreme cases. In the first case, m 2 ^> Am b , the relation ([]]) 
implies that the c quark is nearly on mass shell in the final state X c . Light components of X c 
have an energy of order A and interacting with them the c quark behaves as a heavy quark. In 
particular, for m 2 ^> Am b the momentum px = mxv' of the final state X c can be decomposed 
similar to the momentum p B = m B v of the B meson as 

Px = m c v' + k' , p B = m b v + k, (2) 

with k' and k being the residual momenta. This allows us to apply the heavy quark expansion 
in both m b and m c and introduce the notion of the residual energy Ex = v' ■ k' for the X c state 
similar to that for the B meson, 

£x = m x - m c , 6 = m B - m b . (3) 



In the limit of the "light" c quark, m 2 <C Am;,, there are two different subregions in (H). At m x — 
m 2 c = 0(m c A), or equivalently at ex <C m c , the c quark may be considered as a heavy; we may 
apply HQET and obtain results for an infinitely heavy c quark plus small l/m c corrections. In 
the opposite case, m b A ^> m x ^ m 2 ., we have ex 3> m c , which means that the light components 
of X c have the energy much bigger m c and interacting with them the c quark behaves effectively 
as massless. 

2. Heavy to heavy transitions 

In the region ex <^ m c , the c quark can be treated as heavy and we analyze heavy to heavy meson 
transition B — > WX C in the leading heavy mass limit neglecting both l/m b and l/m c corrections. 
The momenta of the mesons can be decomposed as in (0) and (0) with the residual energy eo 
and ex of order A. We may switch off the spins of the b and c quarks in this approximation, 
and consider the decay b v —> c v *W with b v , c v i and W spinless. Spin Clebsh-Gordan factors can 
be trivially included later. 

2.1. Inclusive structure function 

To simplify consideration we choose the effective interaction lagrangian to be gJW^ + h.c. with 
g the coupling constant and J = c v ,b v . Then the decay matrix elements are A^(b — > cW) = g 
and .M (B — > X C W) = g(X c \J\B). The inclusive differential rate of the decay has the form |J 

dT(B -► X C W) = dT (h -> cW)F(#, e)de (4) 

with e being the energy of the light components of the state X c in the rest frame of the c quark, 

and i? the angle between quark velocities v and v' in (|2]). The structure function is defined in 

HQET as 

F($,e)=J2\(X c \J\B)\ 2 S(e-e x ), (5) 

x c 



where the averaging over B polarizations is assumed and the summation is performed over the 
states X c in HQET with the residual energy ex defined in (J|). The states X c = eg can be 
classified according to the angular momentum and parity, j p , of the light component q. Namely, 
the ground state (S-wave) has j p = \ , and the excited P-wave states have j p = \ and | . 
Exclusive B — > X c decays to S- and P-wave mesons X c are described by the form factors fH| , |7| 



(±~,n|J|±~) = £ (n) (chtf) vfu, 

(\\n\J\\-) = 2r 1 ( ;)(ch^)n' 75 «, 

(l + ,n\J\l-) = V3T$(ch#)v^u, 

where | \ ) is the ground state B = bq meson and where u and u^ are spin | and | meson wave 
functions. Substituting these relations into (|5|) we get the "phenomenological" expression for the 
inclusive structure function in terms of the exclusive form factors |?J 

F(#,e) = §(chtf + 1) ]T ie (n) (ch^)| 2 5(e - e^_) 

n 

+2(chtf - 1) £ \r[f 2 { C ^)\ 2 5{e - s%) (6) 

+ (ch^ + l) 2 (ch^-l)5:|r3 ( ;)(ch^)| 2 ( 5(e-4^ + ) + --- 

n 

where the dots denote the contribution of D- and higher wave states. We follow [ 12| in this simple 



derivation using spinless quarks and W. The structure function F({},e) vanishes identically at 
s < e , where e = £ X L- is the residual energy of the ground state meson. At this point it has 
a 5-peak due to the ground state contribution. Then other narrow peaks follow, separated by 
gaps of order A and corresponding to the excited states in (g). As e grows, the peaks become 
wider, and at e ^> A the structure function F($,e) becomes a smooth function of e determined 
by a hard gluon emission. 

2.2. Relation to Wilson lines 

The structure function (|5D can be represented as a Fourier transformed matrix element of the 
product of two currents, 

F(d,e)= fd 4 xe~ iW - x (B\j\x)J(0)\B) , W = p B - Px • (7) 

The matrix element of the current J = c v ,b v in ([?]) can be simplified in the leading heavy quark 
limit. We recall that for e <^C m c , the b and c quarks behave as heavy particles which move along 
their velocities v and v' , respectively, and interact with light components of the meson. Moreover, 
in the leading power of l/m c and l/m& the only effect of this interaction is the appearance of the 



eikonal phase in the heavy quark wave functions |13| . This phase is given in QCD by a Wilson 
line, P exp(i J c dx ■ A(x)), evaluated along the classical quark trajectory C with the gauge field 
Afj,{x) describing soft gluons [13 j. Then, evaluating the matrix element (B| J*(x) J(0)|B), we 
replace the quark field operators as follows 

6„(0) = $,[0, -oo]a v , b\{x) = e im ^'4$_ t ,[-oo,x] , cv(z)4,(0) = e- im ° v '- x <f> v/ [x,0]D v ,(x) 



where m^ and m c are bare quark masses, a v and a v , are heavy quark creation and annihilation 
operators, and D v > (x) is a free cut propagator of a heavy quark, 

The eikonal phases *$> v [x, — oo] and 3v[x, 0] correspond to the motion the b quark with velocity 
v from — oo to the point x and to the propagation of the c quark from the point to x = v't 
with velocity v', respectively, and are given by 

$ v [x, -oo] = Pexp ( i / dsv ■ A(x + vs) J , $ v >[v't, 0] = Pexp (i / dsv' • A(v's) 

As a consequence, we get the following expression for the matrix element in the leading heavy 
quark limit, 

(B|J t (a;)J(0)|B) = f° i^^-^e-^"™ 6 ^ (B|4 $_„[-oo, ft/] $•[&>', 0]$ w [0, -oo]o„|B) 

J-00 27T 

(8) 
Here, three Wilson lines may be represented as a single Wilson line integrated along the path C 
which goes along v from — oo to point 0, then from point along v' to point tv' and finally from 
point tv' to — oo along — v. Then, for the structure function we get the expression as a Fourier 
transformed Wilson line expectation value 

^e iti£ - £0Ch ^Wc{^,t) (9) 

with the Wilson line defined as 

W c = (B(e )\Pexp(iJ c dx-A(x)^\B(e )) (10) 

Here, in the leading heavy quark limit the state |B(eo)) = a v \B) describes only the light com- 
ponents of the B meson with amputated b quark. We notice that since the structure function 
has a meaning of a cross section and not of a Green function, the gauge fields are ordered in (J|) 
along the path C and not in time. However, on different parts of the path C, the path-ordering 
automatically implies the time-ordering. According to the definition (|9]), the structure function 
F(d, e) is a universal distribution which depends only on the properties of the B meson and not 
on the particular form of the partonic subprocess and the effective interaction lagrangian. The 
^-dependence of the structure function comes, apart from shift of e, from the properties of the 
Wilson line as a functional of the integration path C. 

Relation (H) implies that in order to understand the properties of the structure function out 
of QCD one need to know nonperturbative estimate for the Wilson line expectation value. Let 
us invert (P) and perform the Wick rotation t — > —it, 



X 



W{&, -it) = e fe °( ch,? " 1 ) / ds e- t(£ - £o) F(tf, e) , (11) 

where we took into account that F(i9, e) vanishes for e < e$. Then, using eq. (Bl) we can get 
the phenomenological expression for the Wilson line in terms of the exclusive form factors. Once 
we know the Wilson line expectation value we could extract the structure function and the 
exclusive form factors from (|9]) or (0). For instance, the ground state form factor is related to 
the asymptotic behavior of the Wilson line at large t. Although at present we are not able to 
evaluate a Wilson line nonperturbatively, there is limiting case, the small t limit, in which we 
may effectively use the relation ([TT| ) to extract an information about the transition form factors. 



2.3. Renormalization of the structure function 

Both the exclusive form factors and the inclusive structure function depend on the normalization 
point [i. The \x dependence is introduced by the factorization procedure which one has to per- 
form on the differential rate (|4]) in order to separate contributions of gluons with energy bigger 
and smaller than \x to the perturbative coefficient function, dTo, and to the nonperturbative 
structure function, F, respectively. Each of these functions depends on /x, but this dependence 
is compensated in the physical distribution (§]). 

Relation (|9]) implies that the /i-dependence of the structure function follows from the renor- 
malization properties of the Wilson line. The same is true for the heavy quark exclusive form 
factors, which are defined by the matrix elements of the form (X c \ J(0)|B). In the leading heavy 
quark limit, these matrix elements can be represented as expectation values of Wilson line with 
integration path which goes from — oo to point along the velocity v of the b quark and then 
from to oo along the velocity v' of the c quark. The integration path has a cusp at point 
with the angle t? and, as a consequence, the exclusive form factors obey the RG equation JT5 



f^ + /5(^)^ + r cusp Ki9)V (n) (ch^) = o, (12) 

with the same equation for r^ n \ch'd). Here, r cusp (a s ,i?) is the cusp anomalous dimension known 
up to two-loop order [16|. It is important for us that r cusp (o; s , $) depends on the angle •§ and in 



the small velocity limit it has the following asymptotics flEf : 

r cusp (a s , 0) *=° 2 7cusp (a s ) + 0(tf 4 ) , 7cusp (a s ) = ^C F + (^) 2 C F (c A (^-'Q- N f ^j 

(13) 
Since the integration path C has two cusps at points and tv', the Wilson line Wc entering into 
the definition of the structure function obeys the renormalization group equation [] 

u|- + p(jg)j- + 2r cusp K 0) J w(#, -it) = o . (14) 

We conclude from (H) that the structure function F($, e) obeys the same equation. Then, the 
/i-dependence of the structure function and the form factors, eqs. (|TJ|) and (0), is perfectly 
consistent with the expansion in exclusive channels (§). 

As follows from (|Tl~D, t is conjugated to the energy, e, of light particles emitted to the final 
state X c . Large values of t correspond to small residual energy of the final state in the decay 
B — > WX C and, at the same time, to the "large" size Wilson line ([Hf) . This implies that the 
large t limit is essentially nonperturbative. On the other hand, for small values of t the residual 
momentum of X c is formed by a hard gluon radiation. This suggests that the small t behavior 



of the Wilson line entering (^) and ( |TTD can be analyzed perturbatively. 

Perturbative corrections to the Wilson line in the small t limit do not depend on the states 
involved in its matrix element ([10]), so we may consider its vacuum average for simplicity. Ex- 
panding the path-ordered exponential (|i~0|) in powers of the gauge field and applying the Feynman 



1 It is important to notice that the Wilson line ( |l0[) entering this equation does not involve the time-ordering of 
gauge fields. Had we apply the time-ordering the anomalous dimension should be replaced in (jfj) by r cuS p(a s , i?) + 
r cusp (a s ,?7T -■&). 



rules, we get 
W pert (tf, -H) = 1 + 2a s C F fx 4 - D J ,.^ U) _, 



(2tt 



<J(fc 2 )0(V ' ' 



t> • k v ' ■ k 



l-e 



-t(v'-k) 



0(a : 



(15) 

where the dimensional regularization with D = 4— 2e was used. Integration over gluon momentum 
gives rise to a single pole in 1/e. Performing renormalization of W peTt in the MS scheme, we obtain 



the one-loop expression for the Wilson line |T7] 

fit 



ir, 



pert 






at f v 

2 0& coth •& — 1) log h 2 coth •& / dxx coth a; — 1 — d coth $ 

2 jo 



+ 0{al). (16) 



One easily verifies that this expression satisfies the RG equation (14]). To higher orders of 
perturbation theory, the Wilson line W vcn is a dimensionless function of a s (/j), fit and the cusp 



angle ■&. Then, solving the RG equation (14) we get 

W pert (#,-it) = exp f-2y 2 ^r^p^s^),^) + 7 (a B (2/t),i?) 



;i7) 



where 7(a s , , (9) appears as an integration constant. The expression flT7| ) resumes all leading and 
nonleading logarithmic corrections (a£ log™ t) to the Wilson line W. The one-loop expression for 
j(a s , i?) can be obtained from comparison of ( |i~7|) with (|16D : 



7(«s,^) 



c:v\ 



7T 



Cj7 1 + i? coth § — 2 coth $ / <ia;x coth a; J . 



Although perturbation theory perfectly describes all logarithmic corrections to W, but it fails 
to describe uniquely power corrections in t due to the presence of infrared renormalons |0| in 



perturbative series for the Wilson line [19| 



2.4. IR renormalons in Wilson lines 

The standard analysis of the IR renormalons in physical quantities is performed in the large Nf 
limit with a s A^^ fixed. In this limit only quark loop corrections to the gluon self-energy survive, 
and one assumes that their contribution qualitatively describes the asymptotic behavior of the 
whole series of PT. Let us show that this last assumption is not valid for the structure function. 
The structure function gets large perturbative corrections which are resummed into the Wilson 
line expectation value (0). Perturbative expression (O) implies that the Wilson line exponen- 
tially decreases as the size of the integration path, t, increases. By examining the contribution of 
quark loops to this asymptotic behavior we find that the quark loops contribute to the anoma- 
lous dimensions r cU s P and 7 entering into (|17]) starting at a^Nf order. As a consequence, if we 
represent flT7| ) as exp(— X) then only first two terms, 1 — X, in the expansion of the exponent 
survive in the large Nf limit and their asymptotics has nothing to do with the asymptotics of 
the Wilson line (|i~7l) . Thus the large Nf limit is not applicable for the analysis of the asymptotic 
behavior of the Wilson line. 

Instead of taking the large Nf limit, we perform an "improved" calculation of the one-loop 
Wilson line [19]. First, the argument of the coupling constant in (|T5|) is fixed after one takes into 
account higher order corrections to Wc, to be the transverse momentum of gluons, a s = a B (kt), 



and not the virtuality |20| . Second, we use nonabelian exponentiation theorem [^TJ to rewrite 



the Wilson line as an exponential of the lowest-order result, 

W(-it, 0) = exp (2C F ^ D J j^ 2 5(ky(k )a s (k t ) (^--^J(l- e ~^ 

x exp (tAm(l - ch $)) . (19) 



Here, the first exponent describes the interaction of the eikonal current, created by heavy c and 
b quarks, J e ik(&) = t^t — 7^k)i with soft gluon radiation, while the second exponent is due to 
the heavy quark self-energy correction, E. To the lowest order of PT the self-energy is given by 

„/ ,n ^ a n f d D k a«v 2 , . 

v ; p J (27r) D - 1 k 2 (v-l + vk) v ; 

As was noticed in |22|, |15|, for D = 4 the self-energy S(f •/) contains a linear ultraviolet divergence 
which must be included into the renormalization of the heavy quark masses, mb, c —> rnb,c + Am. 
The last factor in ( |T9"D takes into account the change of the plane wave of the heavy quarks, 
exp(— z(mfcf — m c v') ■ x) with x = —itv', due to renormalization of the heavy quark masses. 
However, before the argument of the coupling constant is fixed, the integral in ( |2(]| ) contains only 
one scale and £(i> ■ Z) vanishes in the dimensional regularization as (I ■ v ) — > leading to Am. = 0. 
That is why self-energy corrections did not appear in the one-loop expression ([15]). 

Let us substitute a s = a s {k t ) = 4- / °° du(k 2 / A 2 )~ u into ( |19D and (|20|) and perform the 
integration over gluon momenta to get 

•" = exp (^0W f du iAtrr( - 2u) /-' dx i^W- + tAm(1 " ch g) ) ^ (21) 

To find Am we evaluate the self-energy correction (p0|) 

E(v ■ I) = A ^^f^ f du T 2 (l - «)r(-l + 2u)(-2v ■ //A) 1 - 2 " (22) 

and notice that the linear divergence of ( PUD corresponds to the singularity of the integrand in 



([22] ) at u = 1/2. Hence, Am is given by p3| 



, Cp r du , , 

Am = -K-^- / , 23 

2(3 J \-2u y } 

where integration is performed at vicinity of u = 1/2. Since (At) 2 " = e -«/ a s(i/t)A) ; ^he exponent 



in ( |2T]) has a form of the Borel transformation. The integration in ( |2i~]) is performed over all 
positive values of the Borel parameters u and the integrand contains the singularities generated 
by the T— function at the points u* = |,1,..., the so-called infrared renormalons. Infrared 
renormalons appear in (|21~D due to singularity of the coupling constant at k 2 = A 2 . At the same 
time, integration in (|2T]) over small values of the Borel parameters, u <^C |, or equivalently over 
large momentum k 2 ^> A 2 gives the result for W vert which after renormalization in the MS scheme 



satisfies the RG equation (|14|) . 

The presence of the IR renormalon in the Wilson line implies, first, that perturbative ex- 
pansion of W is not well defined since we have to specify the prescription for integrating IR 

8 



renormalon singularities. Second, using different prescriptions at u* — \ we get the results for 
the exponent of W which differ in power corrections (At) n with n = 1,2.... This means that 
perturbation theory fails to describe power corrections to the Wilson line and in order to make 
W well defined one has to add nonperturbative power corrections to the exponent of fl2lD , 

W(-it, 0) = W peIt (-it, 0)W nonpeit (-it, 0) . (24) 

Thus, the form of nonperturbative power corrections to W nonper t can be examined by calculating 
the contribution of the IR renormalon to ([H]) • We start with the IR renormalon at u* = \ and 



expand the integrand in (|2T[ ) near u = u* to find its contribution to W as 



Using ( |2lf ) we find that the contribution of the IR renormalon at u* = | to the exponent of 
the Wilson line is compensated by the contribution of the ultraviolet renormalon at u* — | to 
Am. Hence, the leading renormalon contribution to the Wilson line is of order 0(A 2 t 2 ) and 
corresponds to u* = 1. Thus, for the Wilson line (^4j) to be well defined, nonperturbative effects 
should contribute to Wnonpert at the level of 0(A 2 t 2 ) power corrections. Expanding the integrand 
in fl2"T|) near u* = 1 we find the corresponding contribution as 

log W pcrt ~ t ■ + t 2 A 2 sh 2 <? -%- / -*L + G(t 3 ) . (25) 

We can parameterize nonperturbative power corrections to W in terms of hadronic matrix 
elements by performing expansion of the Wilson line W, defined in ([II]), in powers of t: 

W nonpert (-it, ■&) = 1 + ait + -a 2 t 2 + ... , 

where a n = (B\b^ (iv 1 D) n b v \B) . Using the general form of the matrix elements (B\b^D^b v \B) = 
c\V^ and (B\b^D^D u b v \B) = Cig^v + c 3 v^v u and taking into account the equation of motion, 
we get ci = and c 2 = — c 3 = f^l^/3 with /i 2 ,^ = (B\b+D 2 b v \B) being one of the fundamental 
parameters in HQET. As a result, the small t expansion of W is given by 

WWrtHt, #) = 1 + • t + \t 2 ^ sh 2 + 0(t 3 ) . (26) 

Comparing this expression with (p5|) we find an agreement with the IR renormalon analysis. 
Namely, the coefficient in front of t vanishes and the term with t 2 has the same dependence on the 
angle $. Moreover, perturbative 0(t 2 ) corrections to W pert contain infrared renormalon ambiguity 
while nonperturbative 0{t 2 ) corrections to W nonpert have ultraviolet renormalon ambiguity in the 
definition of the matrix element /i 2 ,^. However, the sum of both corrections becomes unique due 
to cancellation of IR and UV renormalons and it defines the nonperturbative parameter /z 2 . 

Being taken separately, W peTt and W noripeTt suffer from the presence of the IR and UV renor- 
malons, respectively. The contributions of IR and UV renormalons cancel each other in the 
product ( pip to make the Wilson line as well as nonperturbative parameters like /i 2 well defined. 
The fact, that the IR renormalons appear in the exponent of PUpert implies that the UV renor- 
malons should exponentiate in PUnonpert as well and using the short distance expansion fl26f) one 
can choose the "minimal" anzatz for nonperturbative part W nonper t{— it, $): 

3 



W QonpeIt (-it, 0) = exp ( -i 2 /4 sh 2 ■& + 0(t 3 ) ) , (27) 



where the expansion in the exponent does not contain a term linear in t. Substituting fl27|) into 
the definition ([§) of the structure function and performing Fourier integration, we obtain the 
nonperturbative anzatz for the structure function 



le-e ch$) 2 \ 2 _ 1 

V / W^l 2^ J' ° -3' 



^nonpert M) = ~^=^ ^Xp ( - ^ ^ '- ) , O 2 = -fj* sh 2 ■& , (21 



which describes a Gaussian distribution of the residual energy e with the width a depending on 
the angle •&. For i? — *■ this function becomes d>(e — £o). This means that the light components of 
the B meson do not interact with outgoing heavy quark as it should be due to the conservation 
of the heavy quark current at zero recoil. 

2.5. Renormalized sum rules for heavy— heavy transitions 

Finally, after substitution of (EJJ) and ( |TTD into (EJj) and ( TTTp we arrive at the sum rule 

POO ( ffl rjh \ 

l o dee-t^Ffre) = exp f -te (ch$ - 1) - 2j -^T cnsp (^a s (k t )) + 7 (0,a B (2/i))J 

x (l + ~nlt 2 sh 2 $ + 0(t 3 )Y (29) 

These sum rule take into account perturbative and nonperturbative effects and generalize the 
Bjorken sum rule 0. With eq.(|6|) taken into account, the l.h.s. of (|29|) is given by infinite sum 
over excited states n. For the lower states n the expansion in @ is well defined and after 
integration in (EjJ) it leads to the power series in t similar to that in (EjJ). But for higher states 
n, corresponding to large residual energy of light particles in the final state X c , the sum in (|6|) is 
determined by hard gluon emission and after integration in ( E9"D it gets logt corrections similar 
to fllTD- The exponential factor in the l.h.s. of (|29|) implements the smooth exponential cut off 
on the residual energy at e — 6q ~ 1/1 This should be compared with a sharp cut off 9(e — 1/t) 
used in 0,§. 

We stress that the sum rules contain two completely independent scales, [i and t, and the both 
sides of ( |29"D have the same dependence on \i and separately on t. The inclusive structure function 
and the exclusive form factors do not depend on t while their dependence on the normalization 
point fi is described by the RG equations ( |T4] ) and (0). We notice that the exponent of the 
integral of the cusp anomalous dimension can be combined with the structure function to produce 
F($, e) at jjl — 2/t. Therefore the sum rule (GDI) can be written in the simpler form 



f" de e-^-^FM, e) = e -*o<**-i)+y(*M2/t)) (i + Iffi sh 2 ^ + 0(t 3 )) . (30) 

Jeo M =2/t V 6 / 

which contains only one scale t. In these sum rules t is a free parameter and by differentiating 
the both sides of (|3~0~1) with respect to t, we can get generalizations of the Voloshin || and the 



"third" 1 10] sum rules. We notice however that the possible values of t in the sum rules are 



l/m c <t<l/A. (31) 



restricted as 



The upper limit for t follows from the condition for the expansion in the r.h.s. of (gbj) and (|30| ) to 
be well defined. To understand the lower limit we recall that the structure function F($, e) was 
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defined before under the additional condition that the c quark should be "heavy" , or equivalently 
the residual energy of the final states X c has to be small, e -C m c . This leads to £ 3> l/m c because 
the residual energy in (p0|) is cut at e — e ~ l/£ with e ~ A. 

Considerable simplifications occur in the sum rules in the small velocity limit, *& — » 0. Let us 
expand the both sides of (|6|) to the first order in d 2 . The ground state form factor is |TT 



£(°)(chtf) = 1 - p 2 (p) (chtf - 1) + . . . 

where p 2 {p) is the slope parameter. The transition form factors to higher S-wave states behave 
as £ ( n ) (ch $) = £ (n ) (l)(chi? — 1) + . . . and may be ignored in (H). Contributions of the transition 
form factors to the P-wave states, r\ii and r 3 / 2 , have the threshold suppression ch^ — 1 in (^), 
and, to 0{d 2 ) order, the form factors may be taken at i? = 0. Moreover, contributions of the 
transitions to D-wave (and higher) states indicated by dots in (|6]) have the threshold suppression 
(ch^ — l) 2 (and stronger), and may be ignored. The anomalous dimension ([H|) behaves in the 
small d limit as 

7 (tf, a.) = tf 2 7o(« s ) + 0{$ 4 ) , 7 o(a B ) = ~C F ^ + 0(a 2 s ) . 

9 71 

Expanding (|29| ) to the $ 2 order, we obtain the generalization of the Bjorken sum rule 
E |r^(l)| 2 e-^+2E |r$(l)|VO = -± V(2/£)+7o(« s (2/£))-^+^ 2 +0(£ 3 ) (32) 

n n 4 Z 

where #1/2(3/2) = £ i/2+(3/2+) ~ £ o are the excitation energies of the P-wave states. Here, the values 
of form factors 71/2(1) and 73/2(1) do not depend on the normalization point /i while the slope 
parameter p is defined at p, = 2/£. 



Differentiating (p2|) in £, we can obtain the Voloshin M and the "third" JTO | sum rules: 



( n )+ « — , . r„\ ^ i„\ s^t 



E l^(l)| 2 ^e"W +2E |r$(l)| a *$e "=' ^ 



^0 - I (p 2 (2/£) + 7o(« s (2/£))) - l&t + (9(£ 2 ), (33) 



E K ( /1(DI 2 («)VO + 2E l^(DI 2 (^e-«' 

n n 

= \&+% (P 2 ( 2 A) + 7o(« s (2/£))) + 0(t) , 
where an arbitrary cut-off parameter £ satisfies condition (]3T|) . 

3. Heavy to light transitions 

Let us consider now the decay of the B meson in the limit of the "light" c quark mass, m 2 -C Ami,. 
As was stressed in sect. 1, for small invariant mass of the final state, m 2 x — m 2 = 0(m c A), the c 
quark behaves effectively as a heavy quark and one can apply the results of the previous section 
to describe the B — > X c + W decay in terms of the structure function of heavy to heavy meson 
transition. As the mass of the final state increases, m^A ^> m? x ^> m 2 , we encounter a new 

11 



situation in which the energy of the light components of X c is much bigger m c and the c quark 
becomes effectively massless. As a consequence, the heavy quark expansion in powers of l/m c is 
not applicable in this case and the decomposition of the momentum px of the final state as well 
as the notion of the residual energy of the final state, ex, do not have a sense. In what follows, 
we neglect the mass of the c quark and analyze the differential rate of the decay B — > X c + W in 
the leading l/m& limit. 

In the rest frame of the B meson the b quark decays into the c quark which has a large energy 
rrib — Wq = 0{mb) and small invariant mass (rribV — W) 2 ~ rn x . The total momentum of the final 
state, px = m B v — W, is equal to the sum of momenta of the c quark and light components of the 
B meson. It is convenient to introduce the light-cone variables p = -y^(p°±p 3 ) and p t = {p l -,p 2 ) 
for an arbitrary momentum p.0 Then, in the rest frame the B meson we have ps = ^(l" 1 ", 1~, 0) 
and W = {W + , W~ , 0). Particles in the final states move close to the light-like direction which 
we choose to be the "— " direction on the light-cone, 

p x = 0(m b ) , p x = 0(A) 

with m x = 2PxPx = 0(rribA). Propagating into the final state, the c quark can produce a jet 
of energetic collinear particles which move in the "— " direction and are accompanied by soft 



radiation including the light components of the B meson ]19[ . This should be compared with the 
case of heavy to heavy transition in which the residual energy e is too small with respect to m c 
for the outgoing heavy quark to create a jet. Thus, in the case of "light" c quark we have to take 
into account a new possibility for the c quark to create energetic jet of particles. This makes the 
analysis more complicated compared to the previous case. 

3.1. IR factorization 

The momentum px of the final state is distributed among collinear and soft particles which 
interact with each other. However since the "— " component of the momentum of collinear 
particles is much bigger the "— " component of the momentum of soft particles, this interaction 
depends only on the "+" components of the corresponding momenta which are of the same order 
for soft and collinear particles. The total momentum of the final state, px, is equal to the sum 
of momentum of the light component of the B meson, EqV, momentum k of soft gluons emitted 
by the b and c quarks in the partonic subprocess and momentum I of collinear particles. Then, 
the phase space for k + component is 

-e v + <k + <p x , 

where the lower limit follows from condition for the total momentum of soft radiation, EoV + k, to 
be time-like, the upper limit corresponds to the extreme case when the total momentum of the 
final state p x = msv + — W + = 0(A) is formed by soft particles while collinear particles have 
zero "+" momentum, l + = 0. 

Performing the analysis of the Feynman diagrams corresponding to the differential rate of the 
decay, we find that the contributions of jet and soft subprocesses to dT{B — ► X c + W) may be 
factorized in the leading l/m& limit into the form |TJ 



d 2 Y rpth + 

def(e)a (p x /v + -e) (34) 



dW 2 dW< 



J -E0 



2 We recall that in the light-cone variables, for arbitrary 4-vectors p and k, we have p-k = p + k +p k + — p • k. 
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where the structure function f(e) takes into account interaction of the b and c quarks with soft 
radiation while ctq{j) x /v + ~~ 6 ) describes interaction of the c quark and collinear particles in the 
final state with momentum l + = p x — ev + . 

The function f(e) is analogous to the structure function of the heavy to heavy transition but 
in contrast with the previous case the velocity v' of the outgoing heavy quark should be replaced 
by "— " light-cone direction in which massless c quark moves into the final state. Then, replacing 
v' by the light-like vector y^ = (0 + ,y~,0) in the definition (^) we find the contribution of the 



soft subprocess to the differential rate as a Fourier transformed Wilson line expectation value [24 



f(e) =v + r %- e^ + y~ W c (v+y-) , (35) 

J-oo All 

with Wc defined in ([H]) for an arbitrary path C. Here, the integration path C goes from — oo to 
point along the b quark velocity v, then along the light-cone "— " direction to point y, and then 
from y to — oo along —v. The structure function f(e) is defined in the limit rrif, — > oo, and its 
properties don't depend on m b . The only place where m b dependence appears is the factorization 
formula fl34]). 

Comparing the definition of the structure functions of heavy-heavy and heavy-light transi- 
tions, eqs.(||) and (|35|), respectively, we find that in both cases the structure function is given by 
a one-dimensional Fourier transformed Wilson line expectation value. The Fourier integration 
is performed along the direction of the outgoing c quark which coincides with a time-like vector 
v' for heavy-heavy transition and with the light-cone "— " direction of the jet for heavy-light 
transition. The argument of the Fourier transformation is equal in both cases to the projection 
of the total momentum of the soft radiation onto the c quark direction. In the heavy-heavy 
transition it is equal to v' ■ (k' — k) with the residual momenta k and k' defined in (§) and (§), 
while in the heavy-light case it is k + = ev + . Since in the heavy-light transition the final state 
X c consists of the jet and the soft radiation, the momentum k is given by k = px — Sqv — I- 

Similar to the structure function F, the function / depends only on the properties of the B 
state and not on particular short distance subprocess. This function describes the distribution 
of the light fields' momenta k + = ev + in the B meson, or equivalently the probability to find 
the b quark with certain light-cone residual momentum projection — k + inside the B meson. 
The factorized expression ([31]) for the differential rate has a simple partonic interpretation with 
f(e) being the distribution function of the b quark inside the B meson and o"o the partonic 
subprocess. In contrast with the distribution function, the partonic subprocess o"o can be calcu- 
lated perturbatively and it gets large perturbative corrections to all order of PT which can be 
resummed |1| . 



3.2. Evolution equation for the structure function 

The heavy quark distribution function has a nonperturbative origin and as a function of the 
factorization scale \i it obeys the GL AP-type evolution equation |23J] . Comparing the definitions 
of the structure functions of heavy-heavy and heavy-light transitions, (§) and fl35|) respectively, 
we notice that there is an essential difference in the definition of the integration path entering 
the Wilson line. Namely, the time-like direction v' in (^|) is replaced by the light-cone "— " 
direction in (^). This suggests to apply the results of the section 2.3 by replacing v'^ by the 
light-like vector y^ = (0 + ,y~, 0). In particular, the angle $ between the vectors v and v' should 
be replaced by the angle between v and y which is equal to infinity for y 2 = 0. However, trying to 
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perform the limit •& — > oo in (|I4"D we notice that the cusp anomalous dimension has the following 



asymptotics |16| in the large "d limit 



$— »oo 



J- cuspids 



ft,. 



gv+ 



7T 



Cf r A S"i2 



N f 



r C usp(« s , i?) = ^r cusp (ft s ) + o(^) 

it \ it / \ y ju ii i 

(36) 
and the RG equation flT4j ) becomes meaningless as $ — * oo. This explains why the structure func- 
tion of the heavy-light transition does not satisfy the RG equation ([14]). For y 2 = the Wilson 
line entering into (EJ3) has additional light-cone singularities which modify the renormalization 



properties of the Wilson line [24, 25, 26 



To find the RG equation for the light-like Wilson line we perform the one-loop calculation 
of Wc using (|T5|) and replacing v' by y. For y 2 = 0, the result of integration contains a double 
pole in e and after renormalization in the MS-scheme we obtain the one-loop expression for the 
light-like Wilson line as |24j] 



Wy 



pert \V "y 



1 + ^C F 

TX 



-L 2 + L- —w 2 

24 



\og[i(/j,v + y~ 



+ 1E 



(37) 



where "— iO" prescription comes from the position of the pole in the free gluon propagator in the 
coordinate representation [p4 |, and je is the Euler constant. It is convinient to absorb the Euler 
constant into the definition of the renormalization parameter /i. Due to additional light-cone 
singularities the Wilson line has a double logarithm of fi to one-loop order. The renormalization 
properties of the light-like Wilson lines have been studied in |^, ESI ^(J . It was shown that the 
Wilson line Wc satisfies the following RG equation 



/ d d \ 

fi— + (3(g)— + r cusp (ft s ) [\og(fxv + y - z0) + log(-^ V + *0)] + T(a s ) J W(v + y- 



with r cusp (ft 



defined in 
r(ft s ) is given by f26| 



r(ft s 



TX 



(38) 
and r(ft s ) some anomalous dimension. Two-loop expression for 



^1 C 



18 108 4 SV V V54 18 



N 



One easily verifies that the one-loop Wilson line (|37[) satisfies the RG equation ([38]). In pertur- 
bation theory the light-like Wilson line is a dimensionless function of fiv + y~ . Solving (|38D we 
find the perturbative expression for W: 



(39) 



- J^^(2T cmp (a 8 (kt))log(k t t) +T(a e (k t ))) +j(a s (l/t))\ . 

The one-loop expression for the integration constant 7(ft s ) can be found from comparison of ([39] 
with (|37T) as 



7(«s 



a s & 2 

Up TX 

tx 24 



0(a 2 ) . 



Expression ( |39| ) should be compared with analogous expression (|T7[) for heavy-heavy transition. 

To find the dependence of the structure function on the renormalization point we use the 

definition (R3f) and perform the Fourier transformation of the both sides of the RG equation 
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with respect to y . Since the anomalous dimension in (^) depends on y , the resulting equation 
for the structure function is the evolution equation [p3, [26|] rather than the RG type equation: 



(4 +/,(s) S /(£)= /-> ,i5(£ - £ ' )/(£,) ' (4o) 

Here, the evolution kernel is given by 

/oo d+ 
— e iet (2r cusp (a s ) logfat - iO) + T(a s )) (41) 

and the "— zO" prescription leads to the spectral property 

P(e) = for e < . 

This is the reason for setting to e the upper limit for e' in ([Up. The kernel P(e) has the 
dimensionality 1/energy, and depends on the dimensional argument e. Therefore it can only 



contain two terms 9(e) je and 6(e). After integration in (j41|) we obtain |24|, |26 



P(e) = -2r cusp (a s ) ( ^ j + ( 2T cnsp (a s ) log - + 7e - r(a s ) J 8(e) , (42) 



where the distribution (9(e)/e) + is defined, for an arbitrary mass scale m, by the formula 




o e Jm e 



for any smooth test function <p(e) which falls quickly enough for the last integral to converge. 
We have (Q(e)/e)™ — (9(e)/e)+ = —S(e)log(m'/m), therefore the kernel (|42j) does not, in 
fact, depend on m. We notice that the evolution kernel contains terms proportional to the 
5(e) function which depend on the renormalization scale fi. In the solution of the evolution 
equation (f40| ) for the distribution function f(e), these terms can be absorbed into e-independent 
factor. Moreover, in the factorized expression ([J4]) for the differential rate the /i dependence of 
the structure function is compensated by the \i dependence of the partonic subprocess a . 

It is important to realize that the renormalization properties of the distribution function ( flUD 
follow from the properties of the Wilson lines. In particular, the definition of the Wilson line 
in ([IT]) does not involve time-ordering of gauge field. Had we start with the same definition ( PB] ) 
of the distribution function but with the time-ordering included, the resulting function would 
satisfy the Brodsky-Lepage evolution equation [2l| for the meson wave function rather than the 



GLAP equation. 

The heavy quark distribution function f(e) defined in (|35D is closely related to the distribu- 
tion function of the b quark inside the B meson, fb/B(x), p4| . 2BJ. Namely, in the leading heavy 



quark limit both function coincide in the end-point region x ~ 1 provided that x = 1 — e/rrib- 
As a consequence, the evolution kernel P(e) coincides with the standard GLAP quark splitting 



function P qq (z = 1 — e/rrib) in the limit z —> 1 |24, ^6|. However, there is an important differ- 



ence between these two distribution functions due to the fact that they have different spectral 
properties. In contrast to the ordinary GLAP equation, the kernel (f42|) does not allow one to 
obtain the anomalous dimensions of the HQET local operators O n = bl(iD + ) n b v . The moments 
of the structure function /~ dee n f(e), after renormalization of f(e), still contain ultraviolet di- 
vergences, because e (having the dimensionality of energy) varies up to infinity. This sharply 
differs from the ordinary QCD case, where integrals in the dimensionless variable x varying from 
to 1 due to spectral properties of /b/B^) and can never contain ultraviolet divergences. 
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3.3. IR renormalons in the light— like Wilson line 

The general solution to the RG equation (^) has a form similar to (24 ) 



W(v + y~) = W pcrt (v + y-)W nonpeTt (v + y-) (43) 

with W / nonpert('W + Z/ _ ) being a nonperturbative boundary value for the Wilson line. To find the 
properties of W nonpcTt (v + y~), we perform the analysis of the IR renormalons in the light-like 
Wilson line. Using (|T9|) and replacing v' by the light-like vector y we obtain 





iAmv ■ y) 



where Am comes from the self-energy corrections to the f-line, eqs. (E^) and (|23"D, while the 
light-cone y-line is protected from these corrections. Integration over gluon momentum leads to 
the expression for the light-like Wilson line 

W^-it) = exp (£^^ r du (At) 3 - ^- 2 "." 1 .-") ) exp (-Am.) . (45) 
\ Po r(l — e) Je sin(7rw) / 

which contains the IR renormalons at u* = ~, 1, . . . . Different choices of the prescription in ( f45| ) 
lead to results which differ in power corrections of the form (At) 2 "* . The leading power corrections 
correspond to the left-most singularity at u* — §. We find however that similarly to the case of 
the heavy-heavy transition the residue of the first exponent of (|45 ) at u — \ is compensated by 
the residue of Am, eq.fl2"3|), at u = | 



log W pwt (-it) = - (^ J -^- u +Am)t+ Oif) = ■ t + Oif) . 

Thus, for u* — \ the contributions of infrared and ultraviolet renormalons to the light-like Wilson 
line cancel each other, so that the leading IR renormalon appears in (|45| ) at u* — 1 

logWWHi) = o • t + t 2 A 2 -%- / -*L + o(t 3 ) . 



We conclude [|19| that resummed perturbation theory generates (but fails to describe uniquely) 
0(A 2 t 2 ) power corrections to the Wilson line Wc- At the same time, this means that nonpertur- 
bative effects should also contribute to 0(A 2 t 2 ) power corrections to make the Wilson line ( f43"D 
well defined. Indeed, performing the small t expansion of the Wilson line we find the expression 
similar to (p6|) 

W nonpeit (-it) = 1 + • t + ^t 2 (46) 

o 

in which the coefficient in front of t vanishes due to the equation of motion for the heavy quark. 

Repeating the arguments in favor of (p7[) we choose the nonperturbative anzatz for the light-like 

Wilson line as 

WnonpertHt) = exp fe* + 0(t 3 )) . (47) 
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Then, after its substitution into (|35|) we find the nonperturbative anzatz for the heavy quark 



distribution function 19 



/non pert (e) = -jL= exp (-^j , a 2 = i/4 , (48) 

which describes a Gaussian distribution around e = with the width a. Since we neglected 0(t 3 ) 
and higher corrections in the exponent (f4"?D, this function does not vanish outside the physical 
region e > —e and it is valid only at the vicinity of e — 0. 

3.4. Renormalized sum rules for heavy— light transitions 

To find the sum rules for the heavy quark distribution function we invert the relation (^5[) and 
express the Wilson line in terms of f{e), 

/oo 
dse~ et f(e). 
-eo 

Now we substitute the expressions fl4~3|), fl39|), and ( f46|) for the Wilson line to get 

^dee^f(e) = exp f -_/ -^(2r cusp (a s (fc t ))log(M) +r(a s (fc t ))) + 7 (a 8 (l/*))J 



Using the evolution equation fl38|) we check that the both sides of the sum rules have the same 
dependence on the renormalization scale //. Putting \x = 1/t we find the renormalized Bjorken 
sum rule 

f" de e~ et f(e) = e 7(Qs(1/t)) [ 1 + ^t 2 + 0(t 3 ) ) (49) 

J-£0 fJ,= l/t V 6 / 

in which the distribution function is defined at the scale \x = 1/t. Here, t is a free parameter which 
implements smooth exponential cutoff on large values of the light-cone momentum component 
k + = ev + . Differentiating QUID in t, we obtain Voloshin and the "third" sum rules. 

4. Conclusion 

In this paper, we have considered properties of the structure functions of inclusive heavy meson 
decays B — > W + X c in two extreme cases of heavy and light c quark corresponding to heavy- 
heavy and heavy-light meson transitions respectively. They are given by Fourier transforms (^) 
and (^) of the diagonal matrix elements of the Wilson line going from — oo to along the heavy 
meson velocity v, then to some point y, and then back to — oo along —v. In the case of heavy- 
heavy transitions, y = v't is time-like, while in the case of heavy-light transitions y is light-like: 
y 2 = 0. As a function of the factorization scale, the former structure function is multiplicatively 
renormalizable (|I4]) , with the anomalous dimension being twice cusp anomalous dimension of 
Wilson lines; the later one obeys the GLAP-type evolution equation fl40|) with the kernel (|42|). 

Several lowest moments of the structure functions are given by simple sum rules, the Bjorken |], 
0, Voloshin || || and "third" |T(J ones. However, being taken in their original form these sum 
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rules turn out to be divergent due to contribution of the states X c with large residual energy. 
Since contribution of these states is essentially perturbative, it should be included to the co- 
efficient functions, and should be therefore excluded from the integrals of the matrix elements 
given by the sum rules. This means that had we perform consistent factorization procedure, the 
integrals over residual energy of the states X c in the sum rules should be cut off in some way at 
the factorization scale /i ]7|, ||, [| |10| . This is also necessary in order to reproduce the correct fi 
dependence of parameters entering the theoretical side of the sum rules. However, no consistent 
method to implement such a cut-off was proposed. 

In this paper, we introduced a soft exponential cut-off on the residual energy of the states 
X c instead of a sharp step-function |7|, ||. We show that with such a cut-off included, the 
sum rules involve the Wilson lines analytically continued to an imaginary separation t or y~. 
At a sufficiently large factorization scale fi, the corresponding Wilson lines can be calculated 
theoretically, as a perturbative expansion plus a series on nonperturbative power corrections. In 
this way, we obtain a consistent renormalized generalization of the Bjorken sum rule for heavy- 
heavy ( |30D and heavy-light (|49|) transitions. The former one considerably simplifies in the small 
$ limit (^). Differentiating this sum rule, we obtain the renormalized generalizations of Voloshin 
and the "third" sum rules fl33l) . 

We have analyzed nonperturbative effects in the structure functions by considering ambiguity 
of perturbative corrections to the Wilson lines to higher orders. The non-abelian exponentiation 
theorem |Tl|] together with the one-loop contribution in which the running coupling constant 



is evaluated at the gluon transverse momentum [20] lead to the improved perturbative expres 



sions (|2l|) and (|45|) for the Wilson line. These expressions contain infrared renormalon poles at 



u* — |, 1. . . on the integration contour over Borel parameter, due to the infrared pole in a s (k t j. 
As a result, perturbative expression for the Wilson line, W peTt , depends on ambiguous integration 
prescriptions for dealing with these singularities. Infrared renormalons induce ambiguities into 
Wp er t at the level of power corrections (At) 2 "*, at which nonperturbative power corrections should 
contribute as well. 

We demonstrate that the uncertainty in the structure functions induced by the leading infrared 
renormalon at u — | in the Wilson line is compensated by the similar uncertainty in the heavy 
quark pole mass. This nicely agrees with the fact that there is no power corrections linear in t. 
The contribution of the next infrared renormalon pole at u = 1 to W pert has the same dependence 
on the cusp angle $ as the nonperturbative 0(A 2 t 2 ) correction to Wnonport proportional to fi 2 . 
The uncertainty induced by this pole must cancel with the ultraviolet renormalon uncertainty in 
the meson matrix element fi^. As a result, the sum of perturbative and nonperturbative (D(t 2 ) 
contributions to the Wilson line becomes unique and it defines the fundamental parameter ft 2 . 

Although other nonperturbative contributions to the Wilson line are possible, the requirement 
of the renormalon ambiguity cancellation leads to the exponentiation of nonperturbative correc- 
tions (^) and (pT7p. If we cut the series in the exponents at the quadratic /z 2 term, we obtain the 
model nonperturbative structure functions fl28|) and (^) describing a gaussian distribution of 
the heavy quark momentum in a heavy meson. These distributions do not vanish exactly in the 
regions required by the support properties. For exact fulfillment of these properties, an infinitely 
many terms in the series in the exponents (p7|), ( [47|) are required. 
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